Dynamics of polymer translocation into an anisotropic confinement 
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Using Langevin dynamics simulations, we investigate the dynamics of a flexible polymer translo- 
cation into a confined area under a driving force through a nanopore. We choose an ellipsoidal 
shape for the confinement and consider the dependence of the asymmetry of the ellipsoid measured 
by the aspect ratio on the translocation time. Compared with an isotropic confinement (sphere), 
an anisotropic confinement (ellipsoid) with the same volume slows down the translocation, and the 
translocation time increases with increasing the aspect ratio of the ellipsoid. We further find that 
it takes different time for polymer translocation into the same ellipsoid through major-axis and 
minor-axis directions, depending on the average density of the whole chain in the ellipsoid, <f>. For <j> 
lower than a critical value <j> c , the translocation through minor axis is faster, and vice versa. These 
complicated behaviors are interpreted by the degree of the confinement and anisotropic confinement 
induced folding of the translocated chain. 
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I. INTRODUCTION 

The transport of biopolymcrs through nanoporcs is a 
challenging problem both for polymer physics and molec- 
ular biology. It has received great attention in both ex- 
perimental and theoretical studies due to its important 
role in many crucial biological processes, such as mRNA 
translocation across nuclear pores, the threading of pro- 
teins through biological membrane channel [j]-{l3j , and 
the viral injection of DNA into host cells [lj-[l7| . Recent 
experiments show that the translocation processes have 
potential biotechnological applications, including rapid 
DNA or RNA sequencing [lElSHl, gene therapy [2l[. 
filtration, and controlled drug delivery [22| . 

Due to the loss of the accessible degrees of freedom, the 
passage of polymer through a nanopore faces a large en- 
tropic barrier and needs an external force to overcome it. 
The driving force can be pro vided by an external applied 
electric field in the pore l23r 37l. a pulling force exerted 

a chemical potential, 
43], longitudinal drag 



on the end of a polymer 
binding particles (chaperones) [3J 

flow in a nanochannel [44| , or geometrical confinement of 
the polymer [iEj . 

Most of previous studies focused on translocation into 
an unconfmed trans side. The dynamics of polymer 
translocation into a confined space through a nanopore 
is of great importance, because it is related to many bi- 
ological processes, such as the bacterial gene swapping 
and DNA prepackaging [46| . Experimentally, Smith et 
al. [47} have found that the rate of packaging is constant 
until about 50% of the DNA is packaged, then dramati- 
cally decreases. In addition, they have also observed the 
pause during the packaging process. Using Langevin dy- 
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namics simulations we have investigated the dynamics of 
polymer translocation into the regions between two par- 
allel plane walls and a fluidic channel [Hj]. Compared 
with an unconfmed environment, the translocation dy- 
namics is greatly changed due to the entropic resisting 
force induced by crowding effect of the partially translo- 
cated monomers. For polymer translocation into a closed 
confinement, experimental results for the rate of packag- 
ing and the pause j47j have been confirmed by stochastic 
rotation dynamics simulations [49j - f5l| and Langevin dy- 
namics simulations [52J. Most recently, using Langevin 
dynamics simulations we have investigated the dynamics 
of polymer translocation into a circular nanocontainer 
through a nanopore under a driving force [53| . We have 
further found that the translocation time distribution 
and the scaling exponent of the translocation time as 
chain length depend on the average density of the whole 
chain in the nanocontainer. Certainly, for the prepack- 
aging of DNA into the virus capsid, electrostatic charges 
should be considered. As noted by Ali et al. [54|, when 
DNA charge is less screened or the Debye length in- 
creases, the packing becomes more difficult and may stop 
midway. 

Interestingly, Ali et al. have already observed that 
capsid geometries play very important role in the pack- 
aging and ejection dynamics of polymers. For a flexible 
polymer a sphere packs more quickly and ejects more 
slowly than an ellipsoid. However, in their simulations 
the lengths of the major axis and minor axis of the ellip- 
soid are fixed, and thus it is still not clear how the as- 
pect ratio (anisotropy) of the ellipsoid affects the packag- 
ing dynamics. In addition, in their simulations polymer 
translocation into the ellipsoid only through its major 
axis is explored. An interesting question is what is the 
difference for polymer packaging into the same ellipsoid 
through its major axis and minor axis? 

To this end, using Langevin dynamics (LD) simula- 
tions, we investigate the dynamics of a flexible polymer 
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FIG. 1: Schematic representation of polymer translocation 
through a nanopore into two-dimensional confined ellipsoid 
under an external driving force. The width of the pore is 
1.6a. The lengths of the semi-axis on x and y directions are 
a and b, respectively 

translocation through nanopores into an anisotropic con- 
finement. We emphasize the influence of the aspect ratio 
of the ellipsoid and translocation direction on the translo- 
cation dynamics. 

The paper is organized as follows. In Sec. II, we briefly 
describe our model and the simulation technique. In Sec. 
Ill, we present our results. Finally, the conclusions and 
discussion are presented in Sec. IV. 

II. MODEL AND METHODS 

In our simulations, the polymer is a coarse-grained 
chain of Lennard- Jones (LJ) particles with the finite ex- 
tension nonlinear elastic (FENE) potential (HJ. The ex- 
cluded volume interaction between beads is generated by 
a short range repulsive LJ potential: 

Ulj = H(-) 12 -(-) 6 }+£, (1) 
r r 

for r < 2 1 / 6 ct and for r > 2 1 / 6 cr. Here, a is the di- 
ameter of a bead and e is the depth of potential. The 
connectivity between neighboring beads is modeled as a 
FENE spring with 

UfeneI?) - -ifcJZg ln(l - r 2 /R 2 ), (2) 

where r is the distance between consecutive beads, k is 
the spring constant and i?o is the maximum allowed sep- 
aration between connected beads. 

We consider a two-dimensional geometry as depicted 
by Fig. [T] where the process of polymer translocation 
into a confinement through a nanopore driven by a trans- 
membrane electric field (or molecular motor for real bio- 
logical system) is shown. The shape of the confinement 



is described by 

x 2 /a 2 +y 2 /b 2 = l, (3) 

where a is the length of semi-axis along the x axis and 
b is the length of semi-axis along the y axis. The con- 
finement is a two-dimensional sphere for a — b and a 
two-dimensional ellipsoid for a / b. 

The setup of the confinement is modeled hard 
shell with a pore of width w = 1.6a, which is formed by 
stationary wall particles within a distance a from each 
other. The entrance permits the translocation of one 
monomer at a time. Between all monomer-wall particle 
pairs, there exits the same short range repulsive LJ inter- 
action as shown above. In Langcvin dynamics simulation, 
each monomer is subjected to conservative, frictional and 
random forces, respectively, 

mh = -V{U L j + Ufene) + F cxt - fv< + Pf , (4) 

where m is the monomer's mass, £ is the friction coef- 
ficient, Vj is the bead's velocity, and Ff- is the random 
force satisfying the fluctuation-dissipation theorem (H3 |. 
The external driving force is expressed as F ext = Fx, 
where F is the external force strength exerted exclusively 
on the beads in the pore, and x is a unit vector in the 
direction along the pore axis. 

In the present work, we use the LJ parameters e, a, 
and bead mass m to fix the system energy, length, and 
mass scales, respectively. The time scale is given by 
tLj = {ma 2 /e) 1 / 2 , and the force scale is e/a. To mimic 
single-stranded DNA, we choose a ^ 1.5 nm [53,1111 and 
the bead mass m 936 amu [59| . Wc set ksT = 1.2s, so 
the interaction strength e corresponds to 3.39 x 10~ 21 J 
at actual temperature 295 K. Thus, this leads to a time 
scale of 32.1 ps and a force scale of 2.3 pN (59[. In our 
simulations, the dimensionless parameters are chosen to 
be Rq = 2, k = 7 and £ = 0.7 unless otherwise stated. 
Corresponding to the range of voltages used in the ex- 
periments @, Hli (3 1 the driving force F is set between 2 
and 15. The Langevin equation is integrated in time by 
the method proposed by Ermak and Buckholz [g3] ■ 

Initially, the first monomer is placed at the pore center, 
while the remaining monomers are undergoing thermal 
collisions described by the Langevin thermostat to reach 
the equilibrium state of the system. Typically, each sim- 
ulation data is the average result of over 1000 successful 
translocation events to minimize statistical errors. Here, 
we defined the successful translocation as the event when 
the whole chain enters into the confinement. 



III. RESULTS AND DISCUSSIONS 
A. Theory 

Although the chain conformation of a polymer confined 
in a slit or a tube has been broadly studied [6ll - [67| . it 
is still a challenging problem for a polymer confined to a 
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closed anisotropic system. For a polymer of chain length 
N confined to an ellipsoid, we define a parameter <5 a ;, to 
characterize the aspect ratio (anisotropy) of the ellipsoid, 

5 a b = b/a. (5) 
The density of beads in the two-dimensional ellipsoid 

is 

= Na 2 /4ab. (6) 

According to the blob picture, a confined chain forms a 
self-avoiding walk consisting of n b blob with size Each 
blob contains g = Cfa/p) D beads, where v-in = 3/4 
being the Flory exponent in two dimensions [63l |64j and 
C is a constant, depending on the temperature and the 
solvent. Then, the density of beads in a blob is 

<t>uob = gv 2 /£,t (7) 

Based on the relation of 4>biob = </>, we obtain the number 
of blobs 

n b = N/g = Ia^V(2^ d -i) = L N <f> 2 . (8) 

Thus, for a polymer with a fixed chain length N, n b is 
determined only by the area of the confinement (nab), 
independent of its shape. Certainly, the above scaling 
arguments can only work as long as the minor axis is 
larger than the blob diameter. 

For polymer confined in an isotropic space, the free 
energy cost is proportional to the number of blobs. In 
units of feflT, it is 

F~N<f> 2 , (9) 

which indicates that the free energy only depends on the 
area of the confinement for a fixed chain length. How- 
ever, for polymer confined in an ellipsoid with the same 
area as the sphere, the blob tends to orient along the 
major axis. Thus, the free energy F should include an- 
other term for the blob orientation which depends on the 
shape of the ellipsoid. Compared with polymer confined 
in sphere, the free energy for polymer confined in ellipsoid 
is higher. For polymer packaging, once the translocated 
chain is confined, the resisting force should be larger for 
the ellipsoid than that for the sphere with the same vol- 
ume. As a result, an ellipsoid packs more slowly than a 
sphere, and with increasing the aspect ratio of the ellip- 
soid the packaging rate should decrease. Certainly, for 
nonequilibrated packaging process the shape dependent 
free energy implies more complicated translocation dy- 
namics. 

During the translocation process, with increasing 
the translocated monomers the entropic resisting force 
f(<fi(t)) increases and depends on the shape of the con- 
finement. Owing to the highly nonequilibrium nature 
of the translocation process and the dependence of the 



shape of the confinement, it is difficult to calculate 
f ((M t)). Here, we can write the translocation time as 

r ~ N a /[F - f(<j>)] = N a /[F(1 - fW/F)], (10) 

where a is the scaling exponent of the translocation time 
r with chain length N, and F(l — f(<f>) / F) is the effective 
driving force. For an unconfined translocation, translo- 
cation time is 

Too ~ N a /F (11) 
for weak driving forces, and then we obtain 

i - — ~ m/F. (12) 

T 

Based on this relationship, we can explore the depen- 
dence of the average resisting force f(4>) on the density 
of the chain <j> and the shape of the ellipsoid. 

B. Numerical results 

In our simulations, 8 a b = 1 corresponds to translo- 
cation into a sphere. The shape of the ellipsoid does 
not change but its orientation changes for the pairs of 
5 ab = 0.125 and 8, S ab = 0.25 and 4, 5 ab = 0.5 and 2, 
respectively. Here, 5 ab > 1 and S ab < 1 denote translo- 
cation into the ellipsoid along the minor-axis and the 
major-axis directions, respectively. 

1. Translocation probability as a function of the driving 
force and the density of the chain 

The translocation probability, Ptrans-, refers to the ra- 
tio of the successful translocation events to all tests at 
given parameters in the simulation. As shown in Fig. 
H{a), with increasing the driving force F, Ptrans increases 
rapidly first, and then slowly approaches saturation at 
larger driving force F. For strong driving forces F > 10, 
the effective driving forces F(l — f ((/>)/ F) are dominant, 
and the shape of the ellipsoid doesn't matter for Ptrans- 
However, for weak driving forces F < 5 the entropic re- 
sisting forces f(<f>) are important, and interestingly we 
observe that Ptrans depends on the shape of the ellip- 
soid. 

To see the shape-dependence of Ptrans more clearly, 
we show the Ptrans as a function of <f> under F = 2 for 
different 8 ab , see Fig. &h). Ptrans is almost constant for 
S a b > 1, independent of <j). However, Ptrans has weak <j> 
dependence and decreases with decreasing 5 a b for S ab < 1. 
Obviously, it is more difficult to translocation through 
major-axis direction of the ellipsoid, implying stronger 
entropic resisting force. This is because at the early stage 
of the translocation process the translocated monomers 
are more confined for 5 ab < 1 compared with the cases 
for S ab > 1. 
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FIG. 2: (a) Translocation probability as function of the driv- 
ing force F for N = 128 and different shapes of the ellip- 
soid 5 a b an d densities of the chain (j>. The solid line with 

— F 

Ptrans = — 1 .8125e 2 .2416 -f 1.0067 is drawn to guide the eyes, 
(b) Translocation probability as function of (f> for N = 128, 
F = 2 and different 5 a b- 
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FIG. 3: Translocation time r as a function of the aspect ratio 
of the ellipsoid 8 a t for N = 128 and F = 2 under the density 
of chain (f> — 0.1 and 0.15. The inset shows the case for the 
strong confinement <j> = 0.42. 



2. Effects of S a b and <f> on the translocation time 

With increasing or decreasing 8 a b from 8 a b = 1 
(sphere), the ellipsoid is more anisotropic. Fig. [3] shows 
the translocation time r as a function of S a b for different 
4>. First, we find that with increasing the anisotropy of 
the ellipsoid, r always increases for different <j>. This 
result indicates that compared with an isotropic con- 
finement (sphere), an anisotropic confinement (ellipsoid) 
with the same volume slows down the translocation ve- 
locity as observed by Ali. et al (49J. 



Figure [3] further demonstrates that it takes different 
time for polymer translocation into the same ellipsoid 
through major-axis and minor-axis directions, depending 
on the average density of the whole chain in the ellipsoid, 
<f>. For tfi = 0.1 and 0.15, translocation into the same el- 
lipsoid along the minor-axis direction is faster. However, 
for <fi = 0.42 the results is opposite as shown in the inset 
of Fig. El 

To give more details, we investigate the translocation 
time as a function of <f> for 8 a b — 0.125, 1 and 8, re- 
spectively. As shown in Fig. @J for all values of (f>, 
translocation into sphere is faster for both F = 2 and 
5. For S a b = 0.125 and 1, r initially increases slowly, 
then increases very rapidly with the increasing <p. How- 
ever, for S a b = 8, t initially increases a little faster com- 
pared with S a b = 0.125. As a main result, we find that 
there exists a critical chain density, denoted as <f> c . For 
4> < 4> c , translocation along the minor-axis direction is 
faster, while translocation along the major-axis direction 
is faster for <j> > <f> c , as shown in Fig. [4ja) for F = 2 and 
Fig. HJb) f° r F — 5- With increasing F, <fi c moves to 
smaller value, such as 4> c w 0.265 for F = 2 and 0.242 for 
F = 5. 

To understand the effect of the anisotropy of the ellip- 
soid on the translocation dynamics, we take into account 
the effective driving forces F(l — f(<p)/F). Based on 
Eq. ([HI), we calculate the average resisting force. As 
shown in the inset of Fig. |4ja) for F = 2, for <fi < <j> c 
the reduced resisting force f(<j>)/F for S a b = 0.125 is 
larger than that for S a b = 8. Thus, the effective driving 
forces F(l — f {<$>)/ F) for 5 a b — 8 is stronger than that for 
S a b = 0.125, leading to a faster translocation for S a b = 8 
as observed. However, for <f> > <fi c the effective driving 
forces for S a b = 0.125 is a little higher, resulting in faster 
translocation for S a b = 0.125. For S a b = 0.125, we find 
that f(4>) ~ ' 45 for (f) < cj) c because the translocated 
monomers don't have enough time to diffuse away from 
the entrance, leading to a high density of the translocated 
chain. For <fi > (j> c it crosses over to /(</>) ~ cj) 1 ' 31 . For 
8 a b = 8, f{<f>) ~ <j) 1A7 is observed. 

In order to shed light on the variation of the aver- 
age resisting force for 8 a b = 0.125 and 8, we give the 
chain conformation right after the translocation for the 
chain N = 128 and different F. As shown in Fig. [5ja), 
for (f> < (j) c the translocated chain is more confined for 
S a b = 0.125 compared with S a b = 8, which leads to a 
stronger resisting force as observed. However, for <f> > <f> c 
the translocated chain forms regular folded chain con- 
formation for S a b — 8 as shown in Fig. EKb), result- 
ing in a stronger resisting force as we calculated. For 
larger driving forces, the translocated chain is easier to 
form folded chain conformation. For cf> = 0.35, there 
even occur quadruple layers in the chain configuration 
for F — 5, compared with triple layers for F = 2. This 
leads to a smaller <p c as shown in Fig. [4] In our previ- 
ous study [48| , we also observed the obvious chain folding 
for polymer translocation into a fiuidic channel through 
a nanopore under the strong driving force. Compared 
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FIG. 4: Translocation time t as a function of the density 
of the chain (j> for N = 128 and different aspect ratio of the 
ellipsoid S a t, under: (a) the driving force F = 2 and (b) the the 
driving force F = 5. The inset shows 1 — (too/0) ~ f(<f>)/F 
as a function of <j>. 
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FIG. 5: (Color online) The chain conformation at the moment 
after the translocation for the chain length N = 128, and 
different density of the chain <j) under the driving force F = 2 
and 5: (a) the aspect ratio of the ellipsoid 5 a b = 0.125 and 
(b) 5 ab = 8. 



with the translocation into a channel with parallel walls 
of the same size as the minor axis of the ellipsoid, the 
translocation into relatively thin ellipsoids is slower and 
is easier to form the folded chain. 



3. Typical translocation process and waiting time 
distributions 

To characterize the typical translocation process, in 
Fig. [6]wc show the plot of the number of beads packaged 
against time for a successful translocation event for F = 
2. For both 6 a b = 0.125 and 8, translocation is faster for 
4> = 0.1 than that for (j> = 0.15. For higher <f>, such as 
4> = 0.35, the translocation velocity is much slower for 
Sab = 0.125 and 8 than that for 6 a b = 1 at the late stage 
of the translocation process. Particularly, it is easier to 
form the folded chain for S a b = 8 compared with S a i, = 1 , 
and it takes longer time for last several monomers to get 
packaged. 

To understand the translocation mechanisms, which is 



considerably affected by the nonequilibrium nature of the 
packaging process, it is necessary to explore the dynamics 
of a single segment passing through the pore into the 
confinement. We numerically calculate the waiting times 
for all monomers in a chain of length N. The waiting 
time of monomer s is defined as the average time between 
the events that monomer s and monomer s + 1 first pass 
through the exit of the pore. 

Fig. [7] shows the waiting time distribution for the 
chain of length N = 128 under different chain densi- 
ties and different 6 a b for F = 2 and 5. For <fr = 0.15, 
the waiting times for s are longer for d a b = 0.125 than 
those for S a i, = 8 for F = 2. This is because the re- 
duced resisting force f(<f>)/F) is stronger for S a b = 0.125, 
and it is important for the effective driving force. How- 
ever, the reduced resisting force f(<p)/F) is almost negli- 
gible compared with the effective driving force for F = 5. 
Thus, there is almost no difference in waiting times for 
S a b = 0.125 and 8. In addition, waiting time distribu- 
tions are almost symmetric with respect to the middle 
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FIG. 6: Number of beads packaged against time for the chain 
length N = 128 and the driving force F — 2. (a) (j> = 0.1 and 
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FIG. 7: Waiting time distribution for the chain length N = 
128 and different </> and S ab . (a) F = 2 and (b)F = 5. 



monomer s — N/2 as previously reported for transloca- 
tion into an unconfined space [23|, [2J, [48| . 

For <f> = 0.35 and F — 2, the waiting times increases 
rapidly after s > N/2 for both 5 ab = 0.125 and 8 due to 
the important role of the resisting forces in translocation 
process. For cf> — 0.42 and F = 5, the waiting times after 
s > N/2 are almost constant for 5 a h = 0.125. However, 
it still has a rapid increase for S a b = 8, resulting from the 
formation of the folded chain conformation. 

These detailed results further confirm the dependence 
of translocation dynamics on the shape of the ellipsoid 
and and the translocation direction. 



IV. CONCLUSIONS 

Using Langcvin dynamics simulations in two dimen- 
sions, we investigate the dynamics of a flexible polymer 
translocation into a confinement under a driving force F 
through a nanopore. Particularly, we address the effects 
of shapes of a confinement on the translocation dynam- 
ics. Compared with an isotropic confinement (sphere), an 
anisotropic confinement (ellipsoid) with the same volume 
slows down the translocation velocity. With increasing 
the aspect ratio of the ellipsoid the translocation time 
increases. We further find that it takes different time 
for polymer translocation into the same ellipsoid through 
major-axis and minor-axis directions, depending on the 
average density of the whole chain in the ellipsoid, (j>. For 
cf> lower than a critical value </> c , the translocation through 
minor axis is faster, and vice versa. These complicated 
behaviors are interpreted by the degree of the confine- 
ment and anisotropic confinement induced folding of the 
translocated chain. 

In the future, it is very necessary to take into account 
the chain stiffness on the translocation dynamics. We 
expect that the anisotropic confinement will play a more 
important role. 
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